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Abstract 

Using the numerical path integral method we investigate the decoherence and relaxation of qubits in spin-boson (SB) and 
spin-intermediate harmonic oscillator (IHO)-bath (SIB) models. The cases that the environment baths with low and medium 
frequencies are investigated. It is shown that the qubits in SB and SIB models have the same decoherence and relaxation as 
the baths with low frequencies. However, the qubits in the two models have different decoherence and relaxation as the baths 
with medium frequencies. The decoherence and relaxation of the qubit in SIB model can be modulated through changing the 
coupling coefficients of the qubit-IHO and IHO-bath and the oscillation frequency of the IHO. 
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I. INTRODUCTION 



A two-level system, namely a qubit in today's mod- 
ern physical language is a most simple quantum physical 
model. If we only consider a closed qubit its dynamics 
is simple enough and has been well known. However, in 
fact, there is not any system being separated from its en- 
vironment, in nature. The interaction between the qubit 
and its environment should be considered when we inves- 
tigate the qubit 's dynamics. The interaction make the 
open qubit model become very complex and it may not 
be solved analytically exactly. In order to investigate the 
model, in most cases, the environment is modeled with a 
thermal bath which is considered being constructed with 
a set of harmonic oscillators. A typical open qubit model 
is the well known spin- boson (SB) model [EHi]. There 
is another open qubit model in which the qubit couples 
to the coordinate X of a harmonic oscillator which we 
shall sometimes call the "intermediate harmonic oscilla- 
tor" (IHO) and which in turn is coupled to a bath 
We call the model the spin-IHO-bath (SIB) model. 

It has be found that the two models have many physi- 
cal correspondences. In recent years the investigations of 
the two models have attached great interest of scientists 
in felids of quantum information, electron transfer, etc. 
0, i, i, 03, [H [H [H El. These are motived by two 
exciting challenges, quantum computation and solar en- 
ergy utilization. It is believed that quantum computers 
may perform some useful tasks more efficiently than their 
classical counterparts. Despite the great promises of per- 
forming quantum computations, however, there are still 
many practical difficulties to be resolved before quantum 
computers might become available in future. One of the 
difficulties is that the qubit has too short decoherence 
time, which is in fact a central impediment for practical 
solid-state qubits to be taken as the cell of quantum com- 
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puters. Many significant results in the field, not only the- 
oretical but also experimental, have been achieved. Most 
of the theoretical research is based on the SB and SIB 
models. On the other hand, the SB and SIB models can 
be used to describe the electron transfer in chemical and 
biological molecules. It has been recent found that the 
coherence is very important to electrons for transferring 
energy in biological systems [l5|, [l6| . So it is important 
to investigate the decoherence and the relaxation of the 
electron spin (qubit). 

In order to investigate the decoherence and relaxation 
of qubit in SB and SIB models, essentiall y, o ne must un- 
derstand the dynamics of the qubit [l7l Il8l [l9| . If the 
qubit energy splitting (denoted by A hereinafter) is not 
equal to zero, the two models are not exactly solvable. 
However, they can be analyzed using adiabatic renor- 
malization in which a systematic weak damping approx- 
imation must be used. They can also be investigated 
with some approximation methods based on the pertur- 
bative scheme which also asks for the system (qubit) 
weakly coupling to its environment. Many other meth- 
ods [20, M, Si SI Si, SI for solving the models have 
been proposed in recent years, most of which are based 
on the Born-Markov approximation. However, it has 
been pointed out that the use of the approximation is 
inappropriate at the large tunneling amplitude and low 
temperatures. So it is important to find out some meth- 
ods to accurately estimate the dynamics of the qubits 
in the two models. Based on the insight into the dy- 
namics we may understand the decoherence and relax- 
ation better and may bring forward some schemes on 
how to suppress them. An excellent, accurate numerical 
method based on the qusiadiabatic propagator path inte- 
gral (QUAPI) [2^, S3 may be a suitable tool for solving 
the two models. To our problems we choose the iterative 
tensor multiplication (ITM) algorithm for the numerical 
scheme. As Makri (2(| S3 addressed that the method is 
non-Markovian and it can make the calculations accurate 
enough even at very low temperatures, large tunneling 
-m&pfflWl&Aimd strong couplings for which the Markovian 
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approximation is unsuitable. In paper I [28| using the 
ITM algorithm we investigated the dynamics and then 
the decoherence and relaxation of the qubit in the SB 
and SIB models as the bath modes with higher frequen- 
cies than the oscillating frequency of the IHO. In this pa- 
per by using the same ITM algorithm we shall investigate 
the decoherence and relaxation of the qubit in low- and 
medium-frequency baths. Where we consider one kind of 
bath, the Ohmic bath interacting with the qubit. The 
low-frequency bath denotes that the Ohmic bath modes 
with lower frequencies than the oscillating frequency of 
the IHO and the medium-frequency bath refer that the 
Ohmic bath modes with the frequencies approaching to 
the IHO's one. 



bath mode is not infinite the effective spectral density is 
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where T = n 2 r)/2M. 
Johm (w) and Jeff (uj) 
plotted in Fig. 1. 



The spectral density functions of 
versus bath modes' frequencies are 



II. MODELS AND DYNAMICS 



The Hamiltonian of the SB model is 



H SB = 7; + Aa x )+J2 



2m,: ~^ 2 ' ' I. 1 Tit;~cj 

(i) 



IB 



+ E 



p 2 



1 



A^) + - + -A^(I 



Act, 



2m, 



—rrijLj' 
2 1 



KCjX 



(3) 



Suppose the bath has an Ohmic spectral density 



J ohm (w) = J E —t (w -un) = l^ue-^ . (2) 
2 i rriiUJi 2 



Here, £ is the dimensionless Kondo parameter [10|, [29[ 
(the relationship of £ with the friction coefficient r\ is 
£ = 2'q/nh (i = x, z) are the Pauli matrix, and 

uj c is the high-frequency cut-off of the bath modes. This 
is a well-known quantum dissipation model and it has 
been widely investigated [l], 0] . 

If we consider the qubit coupling to the coordinate X 
of a single IHO which in turn is coupled to a bath, and 
if we let the couplings be linear, the Hamiltonian of the 
SIB system reads 



Fig. I 

Fig. 1: The spectral density functions J hm ( w ) 
and Jeff (w) versus the frequency U) of the bath 



modes, where A = 5 x 10 Hz, Xk = 1. 



0.01, fi = 10A, T = 0.01 K, 
Hz. 



r = 2.6 x 



i = 

10 11 



Here, we set Xn = 1 and other parameters are same as 
in paper I, namely, A = 5 x 10 9 Hz, £ = 0.01, fi = 10A, 
T = 0.01 K, r = 2.6 x 10 11 Hz. From Fig.l we see that 
the spectral density functions have following character- 
istics. When the frequencies of the bath modes are low 
(0 < uj < 0.1A) the Johm (oj) and J e ff (oj) increase lin- 
early with the bath frequency oj and they are equivalent, 
see Fig. 1. When the frequencies of the bath modes are 
medium (0.1A < uj < 11A) the J hm(w) and J e ff(uj) 
behave in their different ways, see the insetted figures in 
Fig. 1. In particular, in the model SIB, the IHO res- 
onate with some mode of the bath within the frequency 
limits. When the frequencies of the bath modes are high 
(11 A < uj < 100A) the J hm (w) and J e ff (uj) decrease 
with the bath frequency uj. We call the bath modes with 
low- (0 < uj < 0.1A), medium- (0.1A < uj < 11 A) and 
high-frequency (11 A < uj < 100A) the low-, medium- 
and high-frequency baths respectively. The range of fre- 
quencies for the low-, medium- and high-frequency baths 
are plotted in Fig. 2 and they are represented with Low 
F, Medium F and High F in the figure. 



where M and P are the mass and momentum of the IHO, 
and the displacement A characterizes the coupling of the 
qubit to the IHO, and the coupling coefficients of 

the IHO to the bath modes. It is shown that the system 
has a one to one map to the following system [H 



Hsib = - (ea z + Acr^+E 
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According to paper I, as the high-frequency cut-off of the 
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Fig. 2: The sketch map on the low- 
and high-frequency baths. 



medium- 



When the bath modes have high frequencies the dy- 
namics of the qubit in SB and SIB models has been inves- 
tigated in paper I. In this paper we investigate other two 
cases, namely, the cases of low- and medium-frequency 
baths. The IHO is resonance with one of the modes of 
the medium-frequency bath. The length of the memory 
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times of the baths can be estimated by the following bath 
response function 



1 r c 
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(6) 

Here, (3 = 1/ksT where ks is the Boltzmann constant, 
and T is the temperature. It is shown that when the real 
and imaginary parts behave as the delta function 5 (t) 
and its derivative S (t) , the dynamics of the reduced 
density matrix is Markovian. However, if the real and 
imaginary parts are broader than the delta function, the 
dynamics is non-Markovian. The broader the Re[a (t)] 
and Im[a (t)] are, the longer the memory time will be. 
The broader the Re[a (t)) and Im[a (t)] are, the more 
serious the practical dynamics will be distorted by the 
Markovian approximation. The memory time of the ef- 
fective bath is affected by T. The larger the T is, the 
shorter the memory time of the effective bath will be. 
Clearly, the value of the T may be different according to 
the difference of the physical systems. For example, when 
the persistent-current qubit is measured by a dc SQUID, 
the system can be modeled by Eq. (g]) with Eq. ([5]), 
here T = 1/R S C S . Typically^ = 100 ft, C s = 5 pF, 
so T ~ 10 11 Hz, see Ref. j32j. Similar to paper I we 
set r = 2.6 x 10 11 Hz in this paper. In Fig. 3 we plot 
the Re[a (t)] and Im[a (t)] of the low-frequency bath as 
(a) J (uj) = J ohm (^), and (b) J (lu) = J e // (lu) and the 
medium- frequency bath as (c) J (lu) = J h m (u), and (d) 
J (lu) — J e ff (lu) . It is shown that the memory times in 
above four cases are all about r m ss 1.5/ A. 



Fig.3 

Fig. 3: The response functions of the Ohmic bath 
in (a) low and (c) medium frequencies and effec- 
tive bath in (b) low and (d) medium frequencies. 
The parameters are the same as in Fig. 1. The 
cut-off frequencies for the two cases are taken ac- 
cording to Fig. 2. 

The dynamics of the qubit is characterized by the time 
evolution of the reduced density matrix, obtained after 
tracing out the bath degrees of freedom, i.e., 

p (s", s'; t) = Tnatk {s"\ e- mt ' h R (0) j m ' n \s') . (7) 

Thoughout this paper we assume that the interaction 
between system and bath is turned on at < = 0, such 
that the initial density matrix factorizes into its system 
and bath components, and the bath is initially at thermal 
equilibrium (29. [3l|: 



R(0) = p(0)®p bath (0), 



(8) 



where p (0) and p^ath (0) are the initial states of the qubit 
and bath. Here, we calculate the reduced density ma- 
trix p(t) by using the well established ITM algorithm 



derived from the QUAPI. This algorithm is a numeri- 
cally exact one and is successfully tested and adopted in 
various problems of open quantum systems [l(J H3] • 
For details of the scheme, we refer to previous works 
[H, H3. The QUAPI asks for the system Hamilto- 
nian splitting into two parts Hq and H env . Here, we 
take H Q = f (ecr z + Aa x ) and H env = H S b - H , or 
H env = Hsib — Hq. In order to make the calculations 
converge we use the time step At = 0.5/ A, which is 
smaller than the characteristic times of the qubits in the 
systems. 



III. DECOHERENCE AND RELAXATION 

The decoherence is in general produced due to the in- 
teraction of the quantum system with other systems with 
a large number of degrees of freedom, for example the 
devices of the measurement or environment. To measure 
the decoherence one may use the entropy, the first en- 
tropy, and many other measures, such as the maximal de- 
viation norm, etc. (see for example Refs. [33j,l34|). How- 
ever, essentially, the decoherence of an open quantum 
system is reflected through the decays of the off-diagonal 
coherent terms of its reduced density matrix 35]. The 
decoherence time denoted by r-i measures the time of 
the initial coherent terms to their 1/e times, namely, 
p i (n, m) 3 (n, m) = p i (n, m) je. Here, n ^ m, and n, 
m = or 1 for qubits. In the following, we investigate the 
decoherence of the qubit in SB and SIB models via di- 
rectly calculating the evolutions of the off-diagonal coher- 
ent terms, instead of some measure of the decoherence. 
Similar to the decoherence, the relaxation of the qubit 
can also be investigated with calculating the evolutions of 
the diagonal elements of the reduced density matrix. The 
relaxation time is denoted by n, which measures the time 
of an initial state to the final thermal equilibrium state 
through estimating the diagonal terms of the reduced 
density matrix, namely, p i (n, n) — > e~ En @ . In the follow- 
ing calculations we assume that the initial state of the 
environment is p bath (0) = ]J k e^ Mfc / Tr fe (e~^ Mk ) and 
initially the qubit in its maximal coherent state. Here, 
Mk = huikb k bk where b\ (bk) are the create (annihilate) 
operators of the k — th mode for the environment. As 
calculating the off-diagonal element p 12 we let e = 10 A 
which can make the p 12 decay stably. If e — » A the p 12 
will decay with some oscillations, which may affect our 
judgement on decoherence times. The closer the two pa- 
rameters are, the more strongly the matrix elements will 
oscillate. When we calculate p u we choose parameters 
e — A because the oscillations of the p n do not affect 
our judgement on relaxation times from the figures. If 
we choose parameters e = 10 A other than e = A in cal- 
culating the Pn, the relaxation of the qubit cannot be 
shown clearly in the figures in a short time because the 
time n is greatly larger than the time t 2 in our problems. 
So, the reader should note that, the increase of the A and 
e will shorten the decoherence and relaxation times, and 
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the decoherence time Ti and the relaxation time n in 
following figures are not comparable because they are 
plotted in different two sets of the parameters. In the 
ITM scheme, one should at first choose the /c max and as- 
sure that fc max Ai is larger than the effective memory time 
T m of the baths. Fig. 3 and numerical tests tell us that 
the calculations are in fact convergent as fc max ^3. It 
is known that the qubit will show different decoherence 
and relaxation when it has different initial states. In the 
low- and medium-frequency baths, the decoherence and 
relaxation of the qubits in their different initial states are 
different, which is similar to the cases of high-frequency 
bath (see Fig. 3 of paper I), we do not plot the evolutions 
of p 12 and p n in different initial states in this paper. Fig. 
4 plots the evolutions of the p 12 (below) and p n (up) of 
the qubit in low-frequency bath in SB and SIB models. It 
is shown that the qubit in SB and SIB models has almost 
same decoherence and relaxation in this case. 



FigA 

Fig. 4: The evolutions of reduced density matrix 
elements p 12 (below) and p lx (up) in SB and SIB 
models in low- frequency bath. The parameters 
are the same as in Fig. 1. 

Fig. 5 

Fig. 5: The evolutions of reduced density matrix 
elements of p 12 (below) and p lx (up) in SB and 
SIB models in the medium-frequency bath (Xk = 
1, or 1.125). The other parameters are the same 
as in Fig. 1. 

Fig.Q 

Fig. 6: The evolutions of reduced density ma- 
trix elements of p 12 (below) and p xl (up) in SIB 
model in medium-frequency bath in different val- 
ues of f2oi the other parameters are the same as 
in Fig. 1. 

Fig. 5 plots the evolutions of the p 12 (below) and p n 
(up) of the qubit in medium-frequency bath in SB and 
SIB models. It is shown that as kX = 1, the qubit in 
SIB model has longer decoherence and relaxation times 
than the qubit in SB model has. As kX increase to 1.125 
the qubit in SIB model has almost same decoherence and 
relaxation times to the qubit in SB model has, but the 
Pn decays to a bigger equilibrium value. Fig. 6 plots the 
P12 (below) and p xl (up) of the qubit in SIB model in 
different fig of the IHO in medium-frequency bath. It is 
shown that the decoherence and relaxation times of the 
qubit in the SIB model increase with the decrease of the 
oscillation frequency flo of the IHO, which is similar to 
the case that the bath has high frequencies. 



IV. CONCLUSIONS AND DISCUSSIONS 

In this paper we have investigated the decoherence and 
relaxation of a qubit coupled to an Ohmic bath directly 
and via an IHO. In our investigations, we fix the tunnel- 
ing splitting A of the qubit. Two kinds of cases com- 
plementing paper I are discussed. The first is that the 
frequencies of the bath modes are low, and the IHO is far 
off the resonance to the bath modes, i.e. the oscillation 
frequency of the IHO is larger than the high-frequency 
cut-off lu c of the bath modes. The second is that the fre- 
quencies of the bath modes are medium and in the case 
the IHO resonate with some mode of the bath. Here, we 
suppose that, in SIB model, the value of the damping co- 
efficient r of the bath to the IHO is intermediate as sup- 
posed in paper I, which make the bath has shorter mem- 
ory time and it is suitable for our using the ITM based on 
the QUAPI. By using the ITM numerical scheme we cal- 
culated the evolutions of the reduced density matrix ele- 
ments of qubit in the two models. Some new results are 
obtained. (1) When the frequencies of the bath modes 
are low (0 < ui < 0.1 A) the qubit in the SB and SIB 
models has almost the same decoherence and relaxation 
times. Here, we set Xk = 1 other than Xk = 1050 in pa- 
per I. If the damping of the bath to the IHO or (and) the 
IHO to the qubit increase, namely Xk increase, the deco- 
herence and relaxation times of the qubit in SIB model 
will decrease, and vice versa. (2) When the frequencies 
of the bath modes are medium (0.1A < u> < 11 A) and 
Xk = 1, the qubit in the SIB model will have longer deco- 
herence and relaxation times than it has in the SB model. 
As Xk increase to 1.125, the qubit in SIB model will have 
same decoherence and relaxation times to the qubit in 
SB model has, but the p u decays to a bigger equilibrium 
value. (3) The decoherence and relaxation times of the 
qubit in the SIB model increase with the decrease of the 
oscillation frequency fio of the IHO. (4) As point out in 
paper I that the decoherence and relaxation times of the 
qubit in the SB and SIB models will increase with the 
decrease of the e and A, which has also not been plot- 
ted in the paper. The longer decoherence and relaxation 
times are necessary for not only the qubits for making 
the quantum computers but also the electrons for trans- 
ferring energy in biological systems. In order to make 
the qubits or electrons in the SIB model with longer de- 
coherence and relaxation times we may try to make the 
flo and Xk smaller. 

If the qubits in SB and SIB in full-frequency (0 < u> < 
100A), or in low- and medium-frequency (0 < u < 11A), 
or medium- and high-frequency (11A < u> < 100A) 
baths, how about the dynamics of the qubits? This is 
an interesting problem. But, because the effective bath 
in the frequency spectra has so long memory times that 
the ITM algorithm in fact can not be applied for inves- 
tigating the dynamics of the qubit in these cases. So, 
some new methods is expected for the goal. In paper I 
and this paper, we suppose that the qubit couple to the 
Ohmic bath, or the IHO couple to the Ohmic bath. If we 
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set the qubit or the IHO couple to other baths, for ex- 
ample, sub-Ohmic and supcr-Ohmic baths, the method 
used in the two papers is also valid, but in order to make 
the calculations convergence one should recalculate the 
memory times of the baths and choose the k max again, 
so that fc max A£ > r m . 
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